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Abstract 

The Casimir force between a pair of parallell plates filled with ionic particles is considered. We 
use a statistical mechanical approach and consider the classical high temperature limit. In this 
limit the ideal metal result with no transverse electric (TE) zero frequency mode is recovered. This 
result has also been obtained by Jancovici and Samaj earlier. Our derivation differs mainly from 
the latter in the way the Casimir force is evaluated from the correlation function. By our approach 
the result is easily extended to electrolytes more generally. Also we show that when the plates 
are at contact the Casimir force is in accordance with the bulk pressure as follows from the virial 
theorem of classical statistical mechanics. 

PACS numbers: 05.20.Jj, ll.lO.Wk, 12.20.-m, 71.10-w 



Electronic address: johan.hoye@ntnu.no 



I. INTRODUCTION 



It is a pleasure to contribute this work to a festshrift volume for Professor Iver Brevik. 
We have had an extensive collaboration through many years on problems connected to the 
Casimir effect. In our works we have fruitfully utilized methods from different fields of 
research. In particular we have explored the statistical mechnanical aspects of the Casimir 
problem. The present contribution is a work that continues in the statistical mechanical 
direction. 

A pair of metallic or dielectric plates attract each other. This is the well known Casimir 
effect, and it is commonly regarded to be due to fluctuations of the quantum electrodynamic 
field in vacuum. However, H0ye and Brevik considered this in a different way by regarding 
the problem as a statistical mechanical one of interacting fluctuating dipole moments of 
polarizable particles. In this way the Casimir force between a pair of polarizable point 
particles was recovered [l|. To do so the path integral formulation of quantized particle 
systems was utilized [Ij]. Before that this method was fruitfully utilized for a polarizable 
fluid [3]. With this approach the role of the electromagnetic field is to mediate the pair 
interaction between polarizable particles. Later this type of evaluation was generalized to 
a pair of parallell plates, and the well known Lifshitz result was recovered Similar 
evaluations were performed for other situations jl, [6|. 

The statistical mechanical approach opens new perspectives for evaluations of the Casimir 
force. Instead of focusing upon the quantization of the electromagnetic field itself one can 
regard the problem as one of polarizable particles interacting via the electromagnetic field. 
It is found that these two viewpoints are equivalent [l|, 0, 0, 0] ■ 

Metals are materials that have electrons that can be regarded as free. When deriving the 
Lifshitz formula they are regarded as dielectric media that have infinite dielectric constant 
for zero frequency. Jancovici and Samaj realized that it should be possible to evaluate the 
Casimir force for metals by regarding an electron plasma. Thus they consideredparallell 
plates filled with charged particles at low density in a neutralizing background jl, 0, [loj . 
Further they considered the classical case, i.e. the high temperature limit. In this situation 
the Debye-Hiickel theory of electrolytes is fully applicable. Then they use the Ornstein- 
Zernike equation (OZ) equation, and utilize its equivalence with the differential equation 
for the screened Coulomb potential to obtain the pair correlation function. This function is 
used to obtain the local ionic density at the surfaces of the plates. The difference between 
local and bulk densities is attributed to the Casimir force in accordance with the ideal gas 
law. The result obtained coincides with a result for ideal metals in the high temperature 



limit. The latter has been a dispute of controversy The ionic plasma result coincides 
with the one where there is no transverse electric mode at zero frequency. This is also in 
accordance with Maxwell's equations of electromagnetism. 

The ionic plasma has also been extended to the quantum mechanical case by use of the 
path integral formalism from a statistical mechanical viewpoint, and it has been shown that 
magnetic interactions do not contribute in the classical high temperature limit [l^ . 

In the present work we reconsider the ionic plasma in the classical limit. We arrive at the 
the same pair correlation function as in Ref. jSj. But we use a different approach to obtain 
the Casimir force. As we see it, our method better utilizes the methods of classical statistical 
mechanics especially for possible further developments. Thus we use the correlation function 
to directly evaluate the average force between pairs of particles in the two plates and then 
integrate to obtain the total force. This is the method used in Refs. [H, 0]- In this way the 
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result of Ref. [8(] is recovered. A noteable feature of this comparison is that it demonstrates 
that the modification of the density profile at the surface is a perturbing effect that can be 
neglected to leading order by our approach. 

With our approach the evaluations are extended in a straightforward way to electrolytes 
of more arbitrary density. To do so known properties of the direct correlation function is 
utilized. The main change with this extension is that the large distance inverse shielding 
length is modified while the Casimir force remains unchanged for large separations. 

An additional result of our approach is that it is shown that when the plates are at contact 
the Casimir pressure more generally is nothing but the contribution to the bulk pressure 
(with opposite sign) that follows from the virial theorem of classical statistical mechnanics. 



II. GENERAL EXPRESSIONS 



Consider a pair of harmonic oscillators with static polarizability a. They interact via 
a potential 1^3182 where si and S2 are fluctuating polarizations. This interaction creates a 
shift in the free energy of the system. This is easily evaluated to be [H 

1 1 °° 1 

-PF = -- ln[l - [a^r] = -J2 -ic^^r (2.1) 

n=l 

with f3 = l/(kBT) where T is temperature and fc^ is Boltzmanns constant. The last sum 
is the expansion performed in Ref. [J] where the two particles were replaced with two plane 
parallell plates. In the latter case the terms can be interpreted as the sum of graph contri- 
butions due to the mutual interaction ip. The a will represent correlations within each plate 
separately while each ip gives a link between the plates while 2n is the symmetry factor of 
the graphs that form closed rings. With plates the endpoints of each link if) should be inte- 
grated over the plates. In the quantum mechanical case there is also a sum over Matzubara 
frequencies upon which a and ip may depend. 

The parallell plates are separated by a distance a. Due to the interaction there will be 
an attractive force K between the plates. This force is found from [3] 

dF _1 aipa dip 

da (31- {aipy da' ^ ' ' 

The fraction in the middle of this expression represents the graph expansion of the pair 
correlation function with the endpoints in separate plates. These graphs form chains where 
each tp forms a link between the plates. Thus we can write 

K = j /i(r2,ri)^^(r2 -ri)(iri(ir2 (2.3) 

where p is number density, h{r2, ri) is the pair correlation function, and tp'^{r2 — ri) = dtp /da 
with the z-direction normal to the plates. For polarizable particles integral (12.31) will also 
contain integrations with respect to polarizations 

For infinite plates integral (12. 3p diverges, so as usual we will consider the force / per unit 
area which then will be 

/ = J h{k±, Z2, zi)ip'^{k±, Z2 - zi) dk^dkydzidz2 (2.4) 

zi<0,Z2>0 
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where the hat denotes Fourier transform with respect to the x- and y-coordinates. (Here we 
have used J fgdxdy = J f g dk^dky / {2tiY and translational symmetry along the xy-plane.) 
Now we can introduce 

(f' = k\ = k1. + ky, with dk^dky = 2'nqdq. (2-5) 

Further with Z2 = U2 + a and Zi = —Ui we then get 

f = — / h{q, Z2, zi)tlj'^{q, Z2 - Zi)q dqdzidz2 (2.6) 
Ztt J 

ui,U2>0 

An interesting feature of resuh (12.61) or (12. 4p is that it is fully consistent with the virial 
theorem in statistical mechanics. This means that when the plates are at contact for a = 
the Casimir force equals the contribution to the pressure from the virial integral with pair 
interaction tp. With a = translational symmetry is also present in the ^-direction, so we 
have 

/ = p2 J h{\r2-r,\)^'^{\r2-r,\)dxdydz,dz2. (2.7) 

zi<0,Z2>0 

With new variable z = Z2 — zi one can first integrate with respect to Z2 which then will be 
confined to the region < Z2 < z. Thus with dz2 = z we obtain (r = r2 — ri) 

f = P^ j HrWzir) dr = ^J /i(r)rV^(r) dr (2.8) 

z>0 

where symmetry with respect to the x-, y-, and z-directions and with respect to positive and 
negative z is used. (It may be noted that the above is correct if the average of ip is zero. 
Otherwise the pair distribution function 1 + h should be used. But for neutral plates as for 
dielectric plates with dipolar interaction this average will be zero.) 

III. PAIR CORRELATION FUNCTION 

To obtain the correlation function we use the Ornstein-Zernike (OZ) equation 

/i(r2,ri) = c(r2,ri) + j c(r2, r')p(r')/i(r', ri) (ir' (3.1) 
which here has been extended to non-homogeneous fluids. The c(r) is the direct correlation 



function. For week long-range forces [13| or to leading order the c(r) is related to the 
interaction in a simple way 

c(r2,ri) = -/5V^. (3.2) 

For plate separations beyond interparticle distances the ip will be small anyway. For a 
plasma at low density we can write for all r 

c(r2,ri)=c(r) = -/3^, (r = r2 - ri) (3.3) 
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where qc is the ionic charge assuming one component for simphcity. (Here Gaussian units are 
used.) To keep the system neutral a uniform background is assumed. As noted in Ref. ^ the 
OZ-equation is now equivalent to Maxwells equation of electrostatics. The similar situation 
was utilized in Ref. |j| for dipolar interactions. 

Since is the electrostatic potential from a charge one has 

V^cir) = 47r(3ql6{r). (3.4) 

With this Eq. (13. ID can be rewritten as 

- 47r/3g,V(r)$ = -47r5(r - Tq), h{r, Fq) = -pq^^, (3.5) 

where r2 and ri have been replaced by r and ro respectively. In the present case with 
parallell plates the number density is 

( p, z<0 

p(r) = < 0, 0< z<a (3.6) 

p, a < z 

with equal densities p = const, on both plates. By Fourier transform in the x- and y- 
directions Eq. (13. 5p becomes 

-kl-Kl]^ = -An5iz-Zo) (3.7) 




dz'^ / 

where the hat denotes Fourier transform and with = AirjJq^p 

z <0 

kI = { 0, <z <a (3.8) 
a < z. 

The K is the inverse Debye-Hiickel shielding length in the media. Solution of Eq. (13. 7p can 
be written in the form 

+ Bel-', zo<z<0 
$ = 27re^''^° { Ce^'^' + Cle''^ < ;z < a (3.9) 
L)e"«"^, a < z 

where q = k±, q,^ = \/kj_ + k?. (For z < zq the solution is the first line of Eq. (13. 9p where 
the resulting exponent of first exponential has changed sign.) 

With continuous $ and (9$ / dz as conditions, one finds for the coefficient of interest 



J. - 4ge(^--^)" A=(^-^X = f3 10) 

(g, + g)2(l-Ae-2'/«)' W + q) {qu + qY' ^' ' 

With this the pair correlation function for 2:0 < and 2; > a is 

h{k^, z, zo) = -2T^pqlDe-''-'-'~'''\ (3.11) 
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IV. CASIMIR FORCE 



Besides h the ip'^ is needed to obtain the Casimir force /. In accordance with Eq. (13. 3p 
the ionic pair interaction is ip = q^/r. Its full Fourier transform is ip = Anq'^/k'^ which is 
consistent with Eq. (13. 4p . With = k\ + k1 this can be transformed backwards to obtain 

Hk±, z-Zo) = 2'Kql . (4.1) 

q 

This is consistent with solution (13.91) for $. The derivative of (14.11) with respect to z is now 
together with expression (13. lip inserted in Eq. (12. 6p to first obtain {z — zq ^ Z2 — zi = 

U1 + U2 + a) 

f = ^ j {-2'Kj3ql)D{2T^ql) j j e-(«''+'')("^+"^+") c/uiciMs g cig 

00 

00 00 

-idq = -Tr-5 / 1 7-TZ:q'dq. (4.2) 



J {qn + qy 2n/3 J 1- Ae-'^i" 



First one can note that this result is precisely result (3.44) in Ref. This is seen by some 
rearrangement of the latter result with the substitutions Kq ^ k, k ^ q/ k, and d ^ a for 
dimensionality z/ = 3. 

Expression (13.101) and result (14.21 may be simplified further with new variable of integra- 
tion 

q = K, sinh t, dq = k cosh t dt. 

With this we have qn + q = ^(cosht + sinht) = kc*, q^ — q = «e~*, and A = e~^* by which 
the Casimir force becomes 

00 

/ = — - / ^sinh^tcoshtdt (4.3) 



where g(t) = it + 2/tasinht. 

For large separation a only small values of t will contribute, and one can put 

g{t) = (2/ta + A)t and sinh^tcosht = t"^. 

With this and expansion of the denominator the force becomes 

r _ 2C(3) _ knTgS) _ ksTg?,) f 6 \ 

^ 27r/3 (2Ka + 4)3 87ra3(l + 2/(/«a))3 Svra^ ^ /' 

The C(3) is the Riemann C-function, C,{p) = Yl'^=i l/"^^- 

As noted earlier 0, this is the ideal metal result for high temperatures when the 
transverse electric mode is absent. Also one sees that for large a the effective separation 
between the plates is increasesd by twice the Debye shielding length, i. e. a — *• a + 2/k. 
Thus for semiconductors the influence of free ions vanishes due to the increase of effective 
separation for decreasing ionic density. The small conductivity of semiconductors has been 
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an issue of some controversy jl4|. It has been argued that small concentrations of free 



ions in semiconductors should be neglected [15||. However, result (14.11) suggests that lack of 



influence for small ionic concentration is due to increased effective separation for vanishing 



K. 



When the plates are in contact, a = 0, the integral fl3.1ip can be evaluated exlicitly. With 

1 — e~'^* = 4e~^* sinh t cosh t one finds 



oo 



/ = — - / e-2*sinhtrft = (4.5) 







For an ionic system at low density this is precisely the contribution to the pressure (with 
opposite sign) from the ionic interaction (beyond the ideal gas pressure) in accordance with 
the virial integral (12.81) . 



V. ELECTROLYTES IN GENERAL 

For higher densities and lower temperatures the direct correlation function c given by 
Eq. (13. 3p will be modified. However, the crucial point is that for large r ^ oo this expression 
is still valid while for small r there will be changes. On the scale of plate separation this 
change will be a term that can be regarded as a 5-function in r-space such that 

c(r2,ri) = Co(r) + r5(r2 - ri) (5.1) 

where co(r) = —fiq^/r and r is a constant that will depend upon the local density. When the 
local density varies the OZ-equation (13.11) can be regarded as a matrix equation. Multiplying 
it from both left and right with p and adding p on both sides of it the equation after some 
rearrangement becomes 

{l-pc)p{l + hp)=p. (5.2) 
Insertion of expression (15. ip then yields 

(1 - pr - pco)p(l + hp) = p. 

p(l + M = r^^, Pe = -^- (5-3) 

1 - peCo 1 - pr 

Thus the only change in the resulting pair corrrelation function php is that p is replaced by 
an effective density pe on the right hand side. In this way only the inverse shielding length 
is affected by which we get k"^ = AT^jSpeql- -But for large plate separations the Casimir force 
(14. 4p does not depend upon n by which the ideal metal result is generally valid for large 
separations for any electrolyte. 



VI. SUMMARY 

The Casimir force between a pair of parallell plates filled with ionic particles has been 
evaluated in the classical high temperature limit. To do so methods of classical statistical 
mechanics have been used. The pair correlation function is evaluated from which the average 
force between pairs of particles in different plates is found. When the plates are at contact 
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the magnitude of the force equals the contribution to the pressure from the virial theorem. 
This latter result makes the force consistent with bulk pressure. The force found is the same 
as the one found earlier in Ref. [sl for charged particles at low density. There the force 
was evaluated on basis of the difference between surface and bulk densities. By the present 
approach it thus follows that this difference in densities can be neglected to leading order. 
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